Transition from Ferromagnetism to Antiferromagnetism in Gai-^Mn^N 
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Using density functional theory, we study the magnetic stability of the Gai-^Mn^N alloy system. 
We show that unlike Gai-xMn^As, which shows only ferromagnetic (FM) phase, Gai-jMn^N can 
be stable in either FM or antiferromagnetic phases depending on the alloy concentration. The 
magnetic order can also be altered by applying pressure or with charge compensation. A unified 
model is used to explain these behaviors. 



The discovery of ferromagnetism in Mn-containing III- 
V semiconductors has attracted significant attention in 
the last decade because the interesting combination of 
semiconductor electronics and metallic ferromagnetism 
creates great potential for developing functional devices 
that manipulate spin, as well as charge. Among the ma- 
terials that have been investigated, GaMnN is one of the 
most interesting systems. On one hand, original theoret- 
ical predictions suggest that high T c ferromagnetism can 
be achieved in the GaMnN alloy [lj, which has resulted 
in an extensive experimental, as well as theoretical, study 
of its physical properties. On the other hand, available 
experimental data often contradict each other. Some re- 
ports show that high T c ferromagnetism is achievable in 
this system HUSH; 

others show that the magnetic 
coupling of the Mn ions in GaMnN is actually antifer- 
romagnetic (AFM) 6J. The exact nature of the mag- 
netism observed in this system is also still under debate 
IE E3] Some groups suggested that the observed ferro- 
magnetism could be due to secondary phases generated 
during growth, and that GaMnN is a spin glass system 
Other groups argue that no secondary phase is ob- 
served in their single-crystal GaMnN samples that show 
ferromagnetism Q. There are also discussions about 
whether the relevant Mn 3d state has d 4 or d 5 + h charac- 
ter, or whether the d levels are localized in the bandgap 
or inside the valence band [Tof . 

Although many of the issues are still under inten- 
sive study, recent ab initio band structure and total en- 
ergy calculations 0, 0, 0] seem to agree that Mn 3d 
levels are located in the gap, and that the interaction 
between substitutional Mn ions is ferromagnetic (FM) 
at low Mn concentration. Because previous ab initio 
studies also find that pure MnN has an AFM ground 
state 14], an interesting question was raised about how 
the magnetic and electronic properties of Gai-^Mn^N 
evolve as a function of the Mn concentration x. Fur- 
thermore, it is now well known that the ferromagnetism 
observed in Mn-containing GaAs is caused by holes in 
the host valence-band-derived states. It would be im- 
portant to understand how the mechanism changes in 
Mn-containing GaN, where the holes are created in the 
Mn d bands. 

In this paper, using ab initio band structure and to- 
tal energy methods, we study the magnetic properties 
of Gai-^Mn^N in the low and high Mn concentration 



regimes. We also study the effects of pressure and charge 
compensation on the magnetic properties of this system. 
We show that unlike in GaMnAs, where p-d coupling 
induced level splitting at the valence band maximum 
(VBM) is the dominant effect that stabilizes the FM 
phase, direct Mn-Mn d-d coupling mediated by anion p 
states is the dominant effect in GaMnN . We find that 
in GaMnN: (i) The Mn 3d states are inside the band gap, 
in agreement with previous calculations, (ii) At low Mn 
concentrations, the Mn atom tends to have a higher mag- 
netic moment and the magnetic ground state is FM. (iii) 
At higher concentration, the magnetic moment of Mn de- 
creases and the ground state changes to AFM phase, (iv) 
When the material is compressed, the system becomes 
more stable in the AFM phase, even at low concentra- 
tions, (v) When Mn is negatively charged, the system 
turns to AFM, even at low concentrations. We show that 
all these results can be explained by the level crossing of 
Mn d spin-up and spin-down states, which varies with 
Mn concentration and pressure. 

The calculations wer e p erformed using an ab initio 
plane wave basis code [l6|, based on the density func- 
tional theory and using ultrasoft pseudopotentials |17| . 
For the exchange and correlation potential, we used the 
gen eralized gradient approximation of Perdew and Wang 
[l8| . All the structural degrees of freedom (lattice param- 
eters and atomic positions) are optimized by minimizing 
quantum mechanical forces and total energy. The Bril- 
louin zone integration is performed using the Monkhost- 
Pack special k points scheme A large number of 

k points and high cut-off energies for the basis func- 
tions are used to ensure a convergence error within a 
few meV. We considered the zinc-blende alloy and as- 
sumed that the same results also hold for the alloy in 
a wurtzite structure. The disorder of the alloy is taken 
into account explicitly through the special quasirandom 
structure approach |2£}. Supercells up to 64 atoms have 
been used. For AFM calculations involving more than 
two Mn atoms, the sign of the magnetic moments on 
each Mn site are initially distributed randomly. For the 
cases where the AFM configuration was higher in energy 
than the FM phase, several other AFM magnetic con- 
figurations are tested to make sure that the FM phase 
indeed has the lowest total energy. 

We find that unlike Gai-^Mn^As, where the ground 
state is always FM, in Gai_ x Mn x N, the magnetic ground 
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FIG. 1: Schematic model for the hole-mediated stabilization 
of ferromagnetic phase in Gai-iMnjN with low Mn concen- 
tration. In this case, the crystal field splitting is smaller than 
the exchange splitting. The state is the hybridized state 
with majority Mn t2d characeter. For simplicity, the N t^ p 
state is not shown. 



state changes with Mn concentration x. At low Mn con- 
centrations, our calculations show that the FM phase 
of the Mn atoms is more stable, in ag reement with 
other previous theoretical calculations 13']. However, at 
high Mn concentrations, the lowest energy state becomes 
AFM. To understand this interesting behavior, we have 
drawn schematically in Figs. 1 and 2 the Mn d energy 
levels and the coupling between them in the spin-up and 
spin down states. Because the spin-orbit coupling is ne- 
glected in our discussion, only states with the same spin 
and symmetry can couple to each other. We find that 
when Mn atoms are added to GaN, the Mn 3d levels are 
introduced inside the band gap with a d A configuration. 
At each spin channel, due to the crystal field splitting, 
the t2d state is above the state. The crystal field split- 
ting increases when Mn concentration increases due to 
the coupling between the Mn i 2 rf and N t 2p states |21| . 

When the crystal field splitting is relatively smaller 
than the spin exchange splitting, the majority spin states 
will be filled, except for the holes in the tid state, while 
the minority spin states are empty (Fig. 1). In the FM 
configuration, the majority spin state of neighboring Mn 
atoms couple to each other (Fig. la), and the level repul- 
sion pushes one level up and one down. Because the tid 
state is not fully occupied, the level repulsion puts the 
holes in th e higher energy level, and can thus stabilize the 
FM phase |15l l22| . The energy gain in this case depends 
on the position of the Mn-Mn pair and increases with the 
hole concentration at the t^d level. In the AFM config- 
uration the majority spin state of one Mn atom couples 
only to the minority spin states of the other Mn atom 
(Fig. lb) with opposite moment. This level repulsion 



FIG. 2: Schematic model for the stabilization of antiferromag- 
netic phase in Gai-zMn^N with high Mn concentrations. In 
this case, the crystal field splitting is comparable to the ex- 
change splitting. 



push the occupied levels down and the unoccupied lev- 
els up, thus can stabilize the AFM phase. The coupling 
increases when the spin exchange splitting between the 
occupied and unoccupied states decreases and when the 
hole concentration at the tid state decreases. 

At low Mn concentration, because the p-d repulsion is 
weak, the exchange splitting is larger than the crystal 
field splitting (Fig. 1) and the FM interaction is larger 
than the AFM one. This explains why GaMnN has a 
FM ground state at low Mn concentration. When the 
Mn concentration increases, the crystal field splitting in- 
creases due to the larger p-d repulsion. The dispersion of 
the Mild band also increases. When part of the major- 
ity spin tid levels becomes higher than the minority spin 
td state, charge transfer will occur between these two 
states, which leads to a reduced spin exchange splitting, 
as shown in Fig. 2. Because reduced spin exchange split- 
ting will enhance the AFM coupling between the Mn d 
majority state and the minority state, the system will be- 
come increasingly stable in the AFM phase when the Mn 
concentration increases. It is worth noting that when the 
crystal field splitting increases further, only the e<j states 
will be occupied and the system will become paramag- 
netic, with no magnetic moment at the Mn site. 

To test the model discussed above, we plotted in Fig. 
3 the total and the projected density of states (DOS) of 
Gai^Mn.N with x = 0.25 and x = 0.75, which has FM 
and AFM ground states, respectively. From the calcu- 
lated projected density of states (PDOS), we find that at 
x = 0.25, in the FM case, the holes are created in the 
spin-up channel, whereas in the AFM phase, the holes 
are created in both spin channels. The calculated mag- 
netic moment is 3.59 /ib for the FM phase and 3.32 /is 
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FIG. 3: Total (solid) and Mn^ projected (dashed) density 
of states for Gai-jMn^N with different concentrations and 
magnetic configurations. FM is for ferromagnetic and AFM 
for antiferromagnetic. The Fermi energy is at zero energy. 



in the AFM phase. The reason that the FM phase has 
a larger magnetic moment is because the AFM cou pling 
shown in Fig. lb mixes filled and empty d states [23j . 
thus reducing the magnetic moment in the AFM phase. 
At x = 0.75, the increase of the Mn concentration also 
increases the p-d repulsion, leading to a large overlap be- 
tween the majority spin t2d levels and the minority spin 
ed levels. Due to the charge transfer between the major- 
ity spin and minority spin states, the minority ea state is 
partially occupied and the magnetic moment is reduced. 
In the FM and AFM phases, the calculated magnetic mo- 
ments are 1.96 and 2.40 /is, respectively. In this case, the 
FM phase has a smaller moment than the AFM phase, 
opposite to that at low Mn concentration. This is be- 
cause at higher concentration, the charge transfer from 
the spin-up t2d level to the spin-down ed level is larger in 
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the FM phase than in the AFM phase (Fig. 2). We see 
that the calculated results are consistent with our model. 

Our discussion above shows that the change from FM 
to AFM in GaMnN when Mn concentration increases 
is due to the increased crystal field splitting and band 
broadening, which leads to a reduced Mn d-d spin ex- 
change splitting. We notice that the same effect can also 
be simulated by applying pressure (or reducing the lat- 
tice constant). This is because under pressure, the in- 
creased p-d coupling increases the crystal field splitting. 
To test this, we repeated the calculation at x = 0.25, 
but at a lattice constant that is 10% smaller than the 
equilibrium lattice constant. We find that under this 
compression, the system indeed becomes more stable in 
the AFM phase, whereas at its equilibrium lattice con- 
stant, it is more stable in the FM phase. The magnetic 
moment in this case is also reduced, being 2.48 [is for the 
AFM phase and 2.41 hb for the FM phase, similar to the 
case of high concentrations. A similar effect was previ- 
ously observed in the surface of GaN, where the distance 
between Mn atoms is smaller Q. 

The model we discussed above (Fig. 1) also shows that 
the FM phase is more stable when holes are in the Mn d 
bands, whereas the AFM phase will be more stable when 
the holes are filled To test this, we have calculated 
the energy difference AE fm-afm as a function of the 
number of electrons added per Mn atom at x = 0.25. 
The results are plotted in Fig. 4. We find that, indeed, 
the system becomes more stable in the AFM phase when 
the added electron reaches 0.62 per Mn atom. We also 
notice that the lattice constant increases when Mn is neg- 
atively charged because of the larger Coulombic repulsion 
caused by the extra electrons. These results confirmed 
our model. They also suggest that the AFM phase is 
possible in GaMnN if the Mn atoms are compensated 
by donors such as Mn interstitials, N vacancies, or Nc a 
antisite defects. 

In summary, we have explained the magnetic behavior 
of GaMnN using a band structure model. We show that 
because the Mn d states are inside the GaN band gap, 
its behavior is quite different from that of GaMnAs. The 
stability of its magnetic phases can change from FM to 
AFM, depending on the Mn concentration, pressure, and 
charge compensation. 

This work is funded by the U.S. Department of Energy, 
Office of Science, Basic Energy Sciences, under Contract 
No. DE-AC36-99GO10337 to NREL. 
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FIG. 4: Energy difference between ferromagnetic and anti- 
ferromagnetic configurations as a function of the number of 
electrons added per Mn atom in Gao.75Mno.25N. 
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